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CHAPTER  I 


INTRODUCTION 


A  new  feed  configuration  for  microstrip  antennas  has  been  proposed 
by  Pozar  [1].  The  feed  and  a  rectangular  patch  antenna  are  shown  in 
Figure  2.2.  The  feed  consists  of  an  open  ended  microstrip  line  which  is 
located  on  a  dielectric  slab  below  the  ground  plane.  The  microstrip 
antenna  is  formed  on  a  separate  dielectric  slab  above  the  ground  plane 
and  the  two  structures  are  electromagnet ically  coupled  through  an 
electrically  small  aperture  in  the  ground  plane  between  them. 

This  design  is  particularly  advantageous  when  applied  to  millimeter 
wave  monolithic  phased  arrays.  In  this  application  the  associated 
active  elements  such  as  phase  shifters  and  amplifiers  would  be  formed  on 
gallium  arsinide,  which  has  a  high  dielectric  constant  (g^=12.8). 
However,  it  is  preferable  to  mount  the  antenna  elements  on  a  low 
dielectric  constant  substrate  in  order  to  increase  the  bandwidth,  the 
radiation  efficiency  and  the  angle  off  broadside  at  which  scan  blindness 
occurs  [2].  With  the  bi-layered  design,  the  antennas  would  be  located 
on  a  separate  substrate  which  would  yield  optimal  array  performance  and 
eliminate  the  competition  for  surface  space  between  the  antenna  elements 
and  the  active  devices.  In  addition,  the  ground  plane  shields  the 
antenna  half  space  from  spurious  radiation  emitted  by  the  feed  lines  and 
active  devices.  Finally,  aperture  coupling  obviates  problems  associated 


with  probe  feeds  at  millimeter  wave  frequencies,  such  as  complexity  of 


construction  and  large  probe  self-reactances  [1]. 

The  scope  of  this  thesis  is  limited  to  the  analysis  of  a  single 
rectangular  microstrip  antenna  coupled  to  a  microstrip  line  by  a 
rectangular  aperture.  The  goal  of  the  analysis  is  to  accurately  compute 
the  input  impedance  as  a  function  of  frequency,  geometry  and  material 
parameters.  The  analysis  employs  the  moment  method  and  the  resulting 
Sommerfeld  type  integrals  are  numerically  evaluated.  Some  analytical 
results  are  compared  to  experimental  results  to  verify  the  theory. 

Two  approaches  are  pursued  in  the  calculation  of  the  input 
impedance.  The  first  approach  is  to  analyze  the  structure  shown  in 
Figure  2.2  directly  whereas  the  second  approach  involves  extending  the 
feed  line  to  infinity  in  order  to  calculate  the  two-port  S-parameters  of 
the  antenna.  In  the  latter  approach  the  input  impedance  can  be 
calculated  by  simple  transmission  line  theory  when  the  feed  line  is 
terminated  in  an  open  circuit.  Although  the  S-parameter  approach  is 
slightly  less  rigorous  than  the  direct  approach,  since  any  interaction 
between  the  open  circuit  termination  and  the  aperture  is  not  taken  into 
account,  the  effect  of  any  stub  length  on  input  impedance  can  be 
calculated  rapidly  once  the  S-parameters  are  known. 


THEORY 


Formulation  of  the  Problem 


A  schematic  of  the  antenna  and  feed  line  is  shown  in  Figure  2.1 
with  impressed  and  scattered  currents  indicated  on  the  different  antenna 
structures.  In  order  to  simplify  the  analysis  the  ground  plane  and 
dielectric  substrates  extend  to  infinity  in  the  x  -  and  y  -  directions. 
The  electric  surface  currents  on  the  patch  antenna  and  feed  line  and  the 
electric  field  in  the  aperture  are  all  assumed  to  be  y  -  directed  (see 
Figure  2.1a).  By  invoking  the  equivalence  principle  the  aperture  can  be 
closed  off  and  replaced  by  magnetic  surface  currents  just  above  and 
below  the  ground  plane  (see  Figure  2.1b).  The  magnetic  current  above 
the  ground  plane  is  simply  the  negative  of  that  below,  to  ensure 
continuity  of  the  tangential  electric  field  through  the  aperture. 

Denoting  the  space  below  the  ground  plane  (z<0)  as  region  a  and  the 
space  above  the  ground  plane  (z>0)  as  region  b  the  total  electric  and 
magnetic  fields  in  each  region  can  be  broken  into  a  summation  of  fields 
due  to  the  various  currents  as  f'ollows: 


=  E  (J.  )  ♦  E  (J^)  ♦  E  (M  ) 

a  a  me  a  f  as 


=  H  (J.  )  ♦  H  (J^)  t  H  (M  ) 

a  a  me  a  f  as 


(1) 

(2) 
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Figure  2.1  Seheraatie  of  antenna  and  feed  showing  impressed  and 
seattered  currents  for;  (a)  original  problem,  (b)  equivalent  problem 
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=  E.  (j  )  -  E.  (M^) 
b  op  b  3 


=  Hw(3  )  -  H.  (M^) 

b  Dp  b  3 


(3) 

(M 


Each  field  on  the  right  hand  side  of  (1)-(y)  is  the  field  due  to  the 

current  specified  with  the  aperture  shorted.  The  unknown  scattered 

electric  currents  on  the  feed  line  and  antenna  are  J„  and  J 

f  P 

respectively.  The  unknown  magnetic  current  representing  the  aperture 

electric  field  is  M  .  In  region  b  the  negative  sign  on  M  can  be 

3  S 

brought  outside  the  parentheses  because  the  field  is  linearly  related  to 
the  current. 

The  first  step  in  the  moment  method  solution  of  this  problem  is  to 
expand  the  unknown  currents  in  a  finite  set  of  basis  currents  of  unknown 
amplitude.  The  surface  current  density  on  the  antenna  element  is 
expanded  in  a  set  of  piecewise  sinusoidal  (PWS)  basis  functions  which 
vary  in  the  y  -  direction  and  are  constant  in  the  x  direction  as  shown 
in  (5)-(9). 


(8) 
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h  =  L  /(N^+1) 
P  b 


is  the  unknown  coefficient  corresponding  to  the  PWS  basis 
function,  3^.  The  effective  dielectric  constant,  is  taken  from 
Carver  and  Mink  [3].  The  patch  length,  i.e.,  the  resonant  length 
dimension,  and  width  are  given  by  and  respectively.  The  lateral 
offset  of  the  patch  in  the  x-direction  is  x  and  the  antenna  dielectric 
substrate  thickness  is  d^^  (see  Figure  2.2).  The  mode  half  width  is  h^^ 
(see  Figure  2.3  a) . 

Because  the  aperture  is  assumed  electrically  short  the  magnetic 
current  can  be  represented  by  only  a  single  PWS  current  mode.  Noting 
that  the  aperture  is  always  centered  about  the  origin  in  Figure  2.2, 
this  current  can  be  written  as 

M  (x,y)  -  M^P(x,y)  (10) 
s 

where 


1  sin  k  (L  /2  -  Ixl)  ,  -L  /2<  x  <L  /2 

MaP(x  y)  _ _ f£— _ !_1_  V  3p  ap 

^  ’  W  sin  k  L  /2  ’  -W  /2<  y  <W  /2 

ap  ap  ap  ap  ^  ap 

k  <  n/L 
ap  -  ap 


V  is  the  unknown  voltage  across  the  center  of  the  aperture.  The 
o 

aperture  length,  i.e.,  the  long  dimension,  is  and  the  width  is  W^^ 
(see  Figure  2.2).  The  parameter  k^^  must  be  determined  by  comparison  of 
analytical  to  empirical  results,  since  its  value  can  vary  the  aperture 
distribution  from  essentially  trianglar  to  cosinusoidal.  It  will  be 


location  of  additional 


(b) 


first  and  last  test  modes 


(c) 


Figure  2.3  Layout  of  basis  modes  and  those  test  modes  which  do  not 
correspond  to  any  basis  modes  on  the  (a)  antenna  element,  (b)  semi¬ 
infinite  feed  line  terminated  in  an  open,  (c)  infinite  feed  line. 


seen  later  that  a  general  expression  for  k 


can  be  found,  which  is  a 


ap 

function  of  the  material  parameters  and  antenna  geometry  and  yields  good 
agreement  between  measured  and  calculated  results  in  all  cases  studied. 

The  current  on  the  feed  line  is  expanded  in  a  combination  of 
subsectional  and  entire  domain  basis  functions.  The  incident  and 
reflected  currents  are  represented  by  traveling  wave  modes  corresponding 
to  the  fundamental  microstrip  mode  over  the  entire  feed  line.  The  total 
traveling  wave  surface  current  is  then 


*  O’'-'" 


where 


-  •  ,  ,1  - 

■  W/  >'• 

•'rer'”-''*  ■  ^  5/  >'■ 


-W^/2<  X  <Wj./2 
-W^/2<  X  <Wj./2 


The  current  reflection  coefficient,  referenced  to  the  end  of  the  line, 

is  r^.  The  length  of  feed  line  from  the  center  of  the  aperture  to  the 

open  circuited  termination  is  L  and  the  feed  line  width  is  (see 

s  f 

Figure  2.2).  The  term  is  the  exact  propagation  constant  of  the 
fundamental  microstrip  mode.  It  is  calculated  in  a  manner  outlined  by 
Jackson  and  Pozar  [3]. 

_  f 

Writing  J.  in  terms  of  real  functions,  in  order  to  facilitate 
trav 

the  numerical  evaluation  of  subsequent  integrals,  yields 


The  cosine  term  in  (17)  violates  the  boundary  condition  on  the  current 

at  an  open  in  all  cases  except  equal  to  -1.  To  allow  for  general 

solutions  for  r  ,  the  cosine  term  is  truncated  X„/^  (A  denotes  guide 
c  I  r 

wavelengths)  from  the  open  ended  termination.  This  is  shown 


schematically  in  Figure  2.3b.  With  this  condition  imposed  the  incident 
and  reflected  currents  are  altered  near  the  end  of  the  line  and  are 
given  by 


*  jQ(x,y)  -  jJg(x 

.y) 

(18) 

jref(x.y) 

-  r^(J^(x,y)  ♦  jJ 

^(x.y)) 

(19) 

where 

i^cos  k^(y-Lg)y, 

-W-/2<  X  <W./2 

f  f 

-(y  ♦TT/(2k®))<  y  <  L  -TT/(2k®) 

*  e  s  e 

(20) 

l^sin  k^(y-Lg)y. 

-W,/2<  X  <W,/2 

j^(x,y)  . 

f  f 

-y  <  y  <L 

(21  ) 

From  (20)  and  (21),  it  can  be  seen  that  the  traveling  wave  current  mode 
is  finite  in  length.  A  semi-infinite  traveling  wave  current  mode  could 
be  used,  but  this  introduces  an  additional  pole  in  the  integrands  which 
contain  this  mode.  Complexities  associated  with  treating  this  pole  as 


well  as  convergence  problems  are  avoided  by  truncating  the  traveling 
wave  current  mode  in  the  above  manner.  The  semi - inf ini te  line  is 
approximated  by  making  the  traveling  wave  mode  a  few  free  space 


wavelengths  long  and  by  employing  an  appropriate  testing  procedure  to  be 


described  later. 

In  the  vicinity  of  the  aperture  and  open  circuited  termination,  PWS 
basis  functions  are  used  to  represent  non-traveling  wave  currents  which 


arise  due  to  these  discontinuities.  Thus,  the  total  scattered  current 
on  the  feed  line  is 


jf.(x,y)  =  I  I  j^(x,y)  +  r  (j^(x,y)+jJ^(x,y)) 

I  CO  3 

n*  1 


where 


1  Sin 

h^  =  (L^-y^)/(N^>l) 


-W^/2<  X  <W^/2 


is  the  amplitude  coefficient  of  the  n''*^  PWS  current  mode,  J^.  The 
microstrip  feed  width  is  W-.  Parameters  h  and  y  are  respectively  the 

I  3  0 

PWS  mode  half  width  and  the  y  coordinate  of  the  left  edge  the  first 
(n-1)  PWS  basis  mode  (see  Figure  2.3b). 

Substituting  the  expressions  for  J  ,  M  ,  J.  ,  and  J.  given  in 
(5) , ( 1  0) , ( 1  8)  and  (22)  respectively  into  (1)-(4)  and  utilizing  the 
properties  of  linear  operators  yields 


—tot  —  /  "ra 


E  (J^)  -  JE  (5^)  +  I®  E^(j®) 
a  c  a  3  n  a  n 


>  r  [  E  (J^)^jt'  (j^)  ]  ♦ 

c  a  c  a  3  a 


?  rCi^P'i 


mm 
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-  H  (J®)  -  *  r  H  (J®) 

ac  as  nan 

n»1 

(26) 

*  r^C  H^(  3  •  v""  5^(S‘P) 

=tot 

^b 

-  E  (3^)  -  V^P  E  (M^P) 

1-1  ° 

(27) 

qtot 

“b 

-  “b^^iJ  " 

(28) 

The  Individual  field  components  on  the  right  hand  side  of  (25)-(28) 

-b  -ap  -a  ~a 

represent  the  fields  due  to  the  basis  current  indicated  (J^^,  M  ,  J^,  J^, 

and  defined  in  (6) , ( 1 2) , (20) , (21 )  and  (23)  respectively)  in  the 
n 

presence  of  the  grounded  dielectric  slab  only.  To  find  the  unknown 

current  amplitudes  which  occur  in  the  presence  of  the  antenna,  aperture 

and  feed,  three  boundary  conditions  must  be  satisfied.  These  are  zero 

tangential  electric  fields  on  the  patch  and  feed  line  and  continuity  of 

the  tangential  magnetic  field  through  the  aperture.  Continuity  of  the 

tangential  electric  field  through  the  aperture  has  already  been  imposed 

by  the  negative  relationship  between  the  magnetic  currents  above  and 

below  the  ground  plane.  These  three  boundary  conditions  yield  tnree 

equations  to  solve  for  the  three  unknown  currents  J  ,  M  and  J-. 

p  s  I 

Setting  the  tangential  components  of  equal  to  zero  on  the  feed 

line  yields 


{  -  ^  I®  E^j®)  -  r  C  E^M^P) 

‘•.nan  cac  as  a 

n-1 

•  Jon  feed  line 


(29) 


where  the  superscript  t  indicates  tangential  components.  Equating 


tangential  components  of  and  in  the  aperture  gives 

3  O 


-  f  H^j^)  -  r  [  ]  -  V®P[  ] 

nan  cac  as  a  os 

n=1 

N 

>  r  I?  }.  .  (30) 

,,  Ibl  ac  as  in  aperture 

1=1 


“tot 

and  forcing  the  tangential  components  of  to  zero  on  the  antenna 

gives  the  expression  below. 


{-  f  E^j^)  >  E^M^P)  =  0  } 

lfl  1  b  1  b 


on  antenna 


(31 ) 


A  testing  procedure  is  employed  in  order  to  solve  for  the  N, 

i.e.,N  =  N  +  N.  2,  unknown  coefficients  such  that  (29)-(3l)  ane 
a  b 

satisfied.  Each  equation  is  satisfied  in  a  least  mean  square  sense  by 
dotting  it  with  an  appropriate  testing  function  and  integrating  over  the 
structure  where  the  boundary  condition  is  to  be  enforced.  A  Galerkin 
testing  procedure  is  utilized  to  enforce  the  boundary  conditions  at  the 
antenna  and  aperture  to  produce  N  +  1  equations.  The  remaining  N  +  1 
equations  are  obtained  by  testing  with  all  PWS  functions  on  the  feed 
line  in  the  neighborhood  of  the  aperture  and  open  termination.  As  a 
result  there  is  one  more  PWS  test  function  (dashed  PWS  mode  shown  in 
Figure  2.3b)  than  basis  mode  .  This  scheme  enforces  the  boundary 
condition  near  the  terminated  end  of  the  line  only,  however  outside  this 


14 


region  the  scattered  basis  current  assumes  the  form  of  the  fundamental 
quasi-TEM  microstrip  mode  which  automatically  satisfies  the  boundary 
condition  on  the  feed  line.  Furthermore,  by  testing  with  this  procedure 
the  opposite  end  of  the  traveling  wave  mode  is  not  "seen”  by  the  test 
modes  provided  the  traveling  wave  mode  is  terminated  relatively  far  from 
the  aperture  (2-31^).  In  this  way  a  good  approximation  to  the 
infinite  line  is  obtained  without  convergence  problems. 

Using  the  following  notation  to  denote  an  inner  product 


<F.G>3 


//  F*G  ds 
s 


(32) 


and  taking  advantage  of  the  linearity  of  the  inner  product,  (29)-(31) 
upon  testing  become  respectively 


,a  c-y.  ,a 


a  c-y^ 


n-1 


■v=P  .  -O.l.S . 


■h  -  r,[  ]  (3M 

3  3p  D  ap  ^  ,  i  D  i 


1-1 


ap 


1-1 
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The  superscript  x  or  y  on  the  field  quantities  indicate  which  of  the  two 
tangential  components  are  selected  as  a  result  of  dotting  the  field 
against  the  test  function.  The  caret  ( '' )  above  a  current  quantity 
denotes  the  testing  function.  The  letter(3)  cited  after  each  inner 
product  indicates  the  surface  (f=feed,  ap=aperture,  p=patch)  over  which 
the  integration  is  carried  out.  Equations  (33)“(35)  can  be  written  in 
matrix  form  as  follows: 


Cz^][i^]  *  ([z'']+j[z^])r^  -  [T^]v^P  =  -[z'']^j[z®] 

(36) 

ii 
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-Cc^JCi^]  -  (c^+jc^)r^  ♦  (y'*^y^)v^P  +  =  c°-jc® 

(37) 

t 

[Z^JCI^]  >  [T^]V^P  =  0 

(38) 

The  matrix  and  vector  elements  and  scalars  are  defined  as  follows; 
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n  =  1,2,.. .,N, 


1,2 . N 


1  =  1,2 . N, 


It  is  convenient  to  consolidate  some  of  the  quantities  above  into  the 
new  matrices,  vectors,  and  scalars  given  below. 


1  (Z^jz^)] 
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X  N  +1 
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matrix 

(53) 
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N  +1 
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row  vector 
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column  vector 

(55) 
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inc 
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c  -jc 

scalar 

(56) 

[-Z  ♦jz  ] 
m  m 

N  +1 
a 

X  1 

column  vector 

(57) 

The  matrix  represents  the  impedance  matrix  for  an  open  ended 

microstrip  line  by  itself.  Likewise,  [Z^]  is  the  impedance  matrix  for  a 

microstrip  patch  alone.  The  quantities  and  y^  are  the  aperture 

admittances  looking  into  region  a  and  region  b  respectively  of  a  slot  in 

a  ground  plane  with  dielectric  slabs  on  each  side.  The  vectors  [C^ 

tot 

and  [T  ]  are  interaction  vectors  representing  the  coupling  between  the 
feed  line  and  aperture.  Coupling  between  the  antenna  and  aperture  is 


contained  in  [C  ]  and  [T  ].  From  reciprocity  it  can  be  seen  that  the 


elements  of  the  T  matrices  are  related  to  the  C  matrices  by 


b  b 
\  ■  -=i  ■ 


a 

'n 


m  =  n  =  1,2 . ,N, 


(58) 

(59) 


a  t  M 

In  (59)  the  t  element  is  actually  the  (m+1)  element  in  the  vector 
m 


[T  ]  because  of  the  numbering  convention  given  in  (33),  i.e.,  m  runs 


from  0  to  N  ,  whereas  all  other  indices  run  from  1  to  N  or  N.  . 

a  a  b 


Substituting  (53)-(57)  into  (36)-(38)  and  solving  yields 
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1  nc  (>0^  tot  Inc 


(60) 


[if^J  =  *  CT^lv^P) 

tot  tot  inc 

,b 


ri"  ] 

‘  tof* 


-[Z^]“’[T^]V®P 


(61) 

(62) 


a  D 

It  should  be  noted  that  V  equals  -V  ^  from  (11)  and  that  T  is  the 

o  c 


bottom  element  in  [I. 

tot 


Input  Impedance  and  S-Parameter  Analysis 


The  input  impedance  referenced  to  the  center  of  the  aperture  is 
easily  determined  from  the  current  reflection  coefficient,  by 


18 


Z.  =  Z 
in  c 


1  -  r  e 
c 


1  +  r  e 
c 


where  Z^  is  the  characteristic  impedance  of  the  line  and  is  the 


length  of  open  circuited  transmission  line  beyond  the  center  of  the 


aperture. 


It  is  also  of  interest  to  calculate  the  S-parameters  of  the 


aperture  backed  by  the  patch  antenna.  This  can  be  done  with  only  a 


small  modification  to  the  previous  analysis.  The  feed  line  is  extended 


to  infinity  in  both  directions  from  the  aperture.  The  total  traveling 


wave  current  now  becomes  the  sum  of  incident,  reflected  and  transmitted 


wave  components  shown  in  (64). 


and  are  defined  as  before  in  (18)-{21)  where  is  set  equal 


to  zero.  J.  is  given  by 
tran  °  ■' 


^tra,*l’'’y^  =  *  jjgCx.-y)] 


where  J  and  J  are  defined  in  (20)  and  (21)  with  L  equal  to  zero  and 
c  s  s 


T^  is  the  current  transmission  coefficient.  The  layout  of  the  basis 


modes  for  the  S-parameter  analysis  is  shown  in  Figure  2.3c.  The  modes 


are  symmetric  about  the  aperture  center.  The  number  of  unknowns  on  the 
feed  line  now  equals  N  +2.  N  +2  equations  are  obtained  by  testing  with 

3  3 

all  PWS  functions  in  the  vicinity  of  the  aperture.  The  test  PWS 
functions  start  a  half  mode  length  to  the  left  of  the  first  expansion 
mode  and  extend  a  half  mode  length  to  the  right  of  the  last  expansion 
mode  (see  Figure  2.3c).  The  cosine  terms  have  been  truncated  a  quarter 
of  a  wavelength  from  the  aperture  to  avoid  a  current  discontinuity  at 
that  point. 

It  is  evident  from  the  symmetry  of  the  basis  currents  that  the 
matrix  elements  representing  the  coupling  of  the  PWS  test  modes  to  the 
transmitted  traveling  wave  current  are  simply  the  matrix  elements  in 
reverse  order  representing  the  coupling  of  the  PWS  test  modes  to  the 
reflected  traveling  wave  current.  Thus,  very  .'ew  additional 
integrations  are  required  to  carry  out  the  S-parameter  analysis  despite 

r  a  , 

the  introduction  of  an  additional  column  and  row  to  the  matrix 

3  3 

and  an  additional  element  to  the  vectors  [C  ],  [T  ]  and  [V  ]  .  If 

U  V  0  X  1 1  L. 

the  same  number  of  PWS  expansion  currents  are  used  on  the  feed  line  in 
the  two  analyses,  then  only  two  additional  matrix  elements  must  be 
calculated  in  the  S-parameter  analysis  as  compared  to  the  terminated 
line  analysis.  As  will  be  seen  subsequently  most  of  the  additional 
effort  required  for  this  analysis  is  in  the  form  of  matrix  fill 
operations  and  a  couple  of  additional  algebraic  operations. 

Rewriting  (36)-(38)  to  account  for  the  extra  unknown,  T^,  and 
incorporating  the  aforementioned  symmetry  relation  yields 


[Z^]CI^]  ^  ([z‘']+j[z^])r^  ♦  ([z'^l^jCZ®])!^  -  [T^lV^P  (66) 

=  -[z^'l-jCz^] 

-CC^][I^]  -  (c°+jc®)r^  -  (c^-^Jc^)T^  ♦  (y^+y‘^)V^P  (67) 

*  [c"][l"]  =  c^-jc^ 

[Z^^ICI*^]  +  =  0  (68) 

where 


-G  _  C 

”  ^(N  +1 )-m 
a 

(69) 

-s  s 

"  ^(N  -^1  )-m 
a 

(70) 

m=0,1,2,....,N  ,N  +1 

cL  d. 


All  other  quantities  are  defined  as  before  in  (39)-(52)  with  the 
appropriate  change  in  matrix  dimensions  due  to  the  change  in  the  range 
of  index  m. 

As  before,  it  is  worthwhile  to  consolidate  some  of  the  matrices 
above  by  defining 


t4t^ 

r  3 
=  [z 

mn 

1  /  ^  ^  \ 
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1  ,  -C  . 'S 

1  <VJ^m 

)]  N  +2  X 
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N  *2  matrix 
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(71  ) 

^^?ot^ 
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(oJjc^)  1 

1  X  N  +2 
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row  vector 

(72) 

^^?ot^ 

■  t'n  1 

r  It] 

c  '  c 

N  +2  X  1 
a 

column  vector 

(73) 

m  -  0, 1 

f  2 ,  •  •  •  • 

,N  ,N  +1 

With  these  definitions,  (60)-(62)  remain  valid.  It  should  be  noted  from 


(73)  that  r  is  the  next  to  last  element  and  T  is  the  last  element  in 
G  c 

r  I®  1 
"  tot-'" 

Since  the  junction  is  symmetric,  the  S-parameters  referenced  to  the 
center  of  the  aperture  are 


1 1 


(74) 


21 


12 


T 

c 


(75) 


The  input  impedance  can  be  easily  calculated  from  the  S-parameter 
analysis  when  the  feed  is  terminated  in  an  open  circuit.  Examining 
Figure  2.3c  reveals  that  the  relationship  of  the  first  (N  -1)/2  PWS 

Si 

modes  (assuming  N  is  always  odd)  and  the  reflected  traveling  wave 
current  is  identical  to  that  of  the  open  ended  feed  analysis  shown  in 
Figure  2.3b.  Thus,  the  matrix  elements  necessary  to  rigorously 
calculate  the  current  reflection  coefficient  seen  at  the  end  of  an  open 
circuited  microstrip  line  alone  exist  within  and  The 

appropriate  elements  reside  in  the  top  (N  +1)/2  rows  of  [Z  ]  and 

3  u  O  t 

[V  .  ]  and  columns  one  through  (N  -1)/2  and  column  N  -1  of  [Z^  .]. 

inc  3  3  cov 

Letting 


N  =  (N  +1)/2  (76) 

3 

[Z]  =  [z__  (z_+jz_)] 


N  X  N 


(77) 


m  =  0,1,... ,N-1 
n  =  1,2 . N-1 


then  r  is  the  open  circuit  reflection  coefficient  for  a  microstrip 

CO  ^ 

line  alone,  f  is  found  by  solving 

CO 


[Z][I]  -  [V] 


(80) 


With  the  S-parameters  and  the  input  impedance,  referenced  to  the 
aperture  center,  of  an  antenna  terminated  by  an  open  circuited  stub  of 
length  Lg  can  be  easily  found  from 


1  +  r. 


in 


'in 


o  1  -  r.. 


in 


where 


in 


1 1 


^2^21^L 
^  -  ^22^L 
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CO 
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(81  ) 


(82) 

(83) 


Dielectric  Slab  Green's  Functions 


From  (39)-(50)  it  can  be  seen  that  the  fields  due  to  a  given 
current  in  the  presence  of  the  grounded  dielectric  slab  are  required. 


The  fields  can  be  found  by  convolving  the  given  basis  current  agains 


v. 


V 


>  • 
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the  appropriate  Green's  function.  The  required  Green's  functions  are 
the  y  -  directed  electric  field  and  x  -  directed  magnetic  field  due  to  a 
y  -  directed  infinitesimal  electric  dipole  at  (x^,  y^,  d^^)  and  the  same 
field  components  due  to  an  x  -  directed  infinitesimal  magnetic  dipole  at 
(x^,  y^,  0)  for  the  grounded  dielectric  slab  where  the  ground  plane  and 
dielectric  to  air  interface  is  located  at  z=0  and  z=d|^  respectively. 

Figure  2.4  shows  the  nature,  direction,  and  location  of  the  two 
sources  necessary  to  obtain  the  various  Green's  functions.  The  presence 
of  the  dielectric  interface  at  z=d^  gives  rise  to  fields  which  cannot  be 
decomposed  into  a  single  transverse  electric  (TE)  or  transverse  magnetic 
(TM)  mode.  For  the  electric  dipole  of  Figure  2.4a  the  fields  are 

A  A  A  A 

constructed  from  TM  to  y  (A^y)  and  TM  to  z  (A^z)  vector  potentials  and 
analogously  the  fields  due  to  the  magnetic  dipole  of  Figure  2.4b  are 

A  A  A  A 

formed  from  TE  to  x  and  TE  to  z  (F^z)  vector  potentials. 

The  vector  potentials  are  derived  using  the  Fourier  transform 
method  which  is  suitable  for  structures  unbounded  in  one  or  more 
dimensions.  By  Fourier  transforming  the  scalar  wave  equation  for  each 
of  the  above  vector  potential  components  with  respect  to  x  and  y,  the 
problem  is  reduced  to  a  one  dimensional  Green's  function  problem  in  z. 
The  X  and  y  dependency  of  the  fields  are  obtained  by  performing  the 
inverse  Fourier  transform. 

In  the  case  of  the  y  -  directed  electric  dipole  of  Figure  2.4a  the 
scalar  differential  equations  which  must  be  solved  are 
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Figure  2.4  Orientation,  location  and  type  of  current  sources  used  to 
solve  for  the  Green's  functions  in  region  b.  (a)  J-source  (b)  M-source 
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Figure  2.5  Orientation,  location  and  type  of  current  source  used  to 
solve  for  the  Green's  functions  in  region  a.  (a)  J-source  (b)  M-source 
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[7^  +  k^]A^2^^*y'^)  “  2;  >d^  (84d) 

where 

k^  »  (85) 

The  free  apace  permittivity  and  permeability  are  and  respectively. 
The  complex  dielectric  constant,  which  takes  into  account  dielectric 
loss,  is  e^.  Fourier  transforming  (8Ua)-(84d)  with  respect  to  x  and  y 
gives 
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Im{k2}  <0,  Re{k2)  >  0 


.2  ,2 

k  +  k 
X  y 


(90) 
(91  ) 


The  quantities  k  and  k  are  the  two  independent  separation  constants  in 
X  y 


the  separation  equations  given  in  (89)  and  (90).  The  sign  on  the 


imaginary  part  of  k^^  and  k^  is  chosen  to  guarantee  convergence  of  the 


integrals  in  (87)  and  (88)  and  the  sign  on  the  real  part  is  necessary  to 
remain  on  the  same  branch  as  determined  by  the  sign  on  the  imaginary 
part  [5].  Solving  (86a)-(86d)  subject  to  the  boundary  conditions  yields 
the  magnetic  vector  potential  components  in  the  spectral  domain. 
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The  transformed  electric  and  magnetic  field  quantities  of  interest  are 
easily  obtained  from 
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Thus,  the  required  Green’s  functions  in  the  space  domain  due  to  the  y  - 

directed  infinitesimal  dipole  are:  (1)  the  y  -  directed  electric  field 

at  z-d^ 
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1  D  D 


1/2 

where Z  =(u  /e  )  and  (2)  the  x  -  directed  magnetic  field  at  z-0 
o  o  o 
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The  above  Green's  functions  are  appropriate  for  finding  fields  in 


region  b  due  to  the  currents  on  the  antenna  element.  The  Green's 


functions  necessary  to  find  the  fields  in  region  a  due  to  the  feed  line 


currents  are  determined  from  the  configuration  of  Figure  2.5a.  The 


fields  in  this  case  are  easily  determined  from  (100) -(  103)  by  a  simple 


b  a 

coordinate  transformation  and  replacing  d.  with  d  and  e  with  e  . 
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Thus,  the  y  -  directed  electric  field  at  z»-d  is 
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go 
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Qp,  (k  ,k  )  «  Qp  (k  ,k  )!  .  (see  eqn.  101) 
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and  the  x  -  directed  magnetic  field  at  z-0  is 
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Qu,  (k  ,k  )  -  -Qu,  (k  ,k  )|  .  (see  eqn.  103) 
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In  the  case  of  the  x  -  directed  magnetic  dipole  at  z-0  in  Figure 


2.4b  the  same  procedure  as  above  is  followed  to  solve  for  the  electric 


vector  potential.  The  Fourier  transformed  differential  equations  for 


the  electric  vector  potential  components  are 


d^  2  -  ■'*^‘^x’'o''‘^vyo^ 

dz  y 


^^2  ^  ‘^Ib^'^zl  (l^x’l^y'^)  “  0* 

0<  z  <d, 
b 

(108b) 

2 

dz  ^ 

z  >d. 
b 

(108c) 

2 

[  2  "■  1^2^^z2(‘^x’1^y’^)  “  0* 

dz  ^ 

z  >d. 
b 

(108d) 

where 

Solving  (1Q8a)-(108<l)  and  applying  the  boundary  conditions  gives  the 
electric  vector  potential  components  in  the  spectral  domain  given  below. 


I  - — -  cos(x,^z)  nil) 

1  b  m 


sin(k  z)  -j(k  x  +  k  y  ) 

]  e  r° 
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F  (k  ,k  ,z) 
z1  X  y 


-i)  -j(k  X  +  k  y  ) 

xr  ,  ^  xo  v'o 

— r — r —  sin(k,.z)  e  ^ 

^b  .j,b  1b 
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F  (k  ,k  ,z) 
x2  X  y 
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m 


(113) 


r.jlK^.ky.x) 


"x’V’^l''(Sb'‘b>  -J'<2<^-‘>b>  -J''<x’<o  *  , 
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The  two  field  quantities  desired  are  obtained  from 


E  (k  ,k  ,z) 
y  X  y 

H^(k  ,k  ,z) 
X  X  y 


_  > 
dz  X 


Jk  F 
X  z 


T-!-  C(e\^-k^)F 
r  o  X  X 


(116) 


From  ( 1 1 1 )  ,  ( 1 1 2) , (1 1 5)  and  (1l6)  the  necessary  Green's  functions  in  the 

A 

space  domain  due  to  the  x  -  directed  infinitesimal  magnetic  dipole  can 

A 

be  readily  found.  These  are:  (1)  the  x  -  directed  magnetic  field  at  z 

-  0 


dk  dk  (117) 
X  y 


“HHxx‘'<x'‘‘y>  ■  Z? 


b,  2  .2, 


4u"k  Z  k,  J  T 
00  1b  e  m 


(118) 


X  {k^^k2(e^>1)sin(k,^d^)co3(k^^djj)  ^  j ( G\^3in^(k^ ^d^)-k;^^^co3‘'(k^ ^d^^) ) }  ] 


(2)  the  y  -  directed  electric  field  at  z  *  d 


j(k  (x-x  )+k  (y-y  )) 

*  X  °  dk  dk  (119) 

X  y 


qEw  (k  ,k  )  =  -Q^,  (k  ,k  ) 
EMyx  x  y  HJxy  x’  y 


(see  eqn.  103) 


(120) 


Equation  (120)  follows  from  reciprocity. 

As  before,  the  Green's  functions  appropriate  for  region  a  shown  in 
Figure  2.5b  can  be  formed  from  those  given  in  (117)“(120).  These 
Green's  functions  are  the  x  -  directed  magnetic  field  at  z  «  0 
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O  O 


II 


j(^jx-x)*kjy-y)) 


dk  dk  (121) 
X  y 


^HMxx^^*^^  =  ^HMxx^^'^^l  b  a 

V  "  V*  %  " 


(122) 


and  the  y  -  directed  electric  field  at  z  =•  -d 


a  if  ,  j('<  (x-x  )+k  (y-y  )) 


dk  dk  (123) 
X  y 


^EMyx^^*^^  °  ■'^Lyx^^'‘'y^l  b  a 

c  -►  e  .  d^  -*•  d 
r  r  b  a 


(see  eqn.  120)  (124) 


Matrix  Element  Expressions 


With  these  Green's  Functions  the  matrix  elements  in  (39)-(50)  can 

be  fully  defined.  In  (39)  E^(J^)  is  the  y  -  component  of  the  electric 

field  at  z  =  -d^  due  to  the  basis  current  in  the  presence  of  the 

dielectric  slab.  This  field  is  obtained  by  convolving  the  Green's 

function  given  in  (104)  against  as  shown  below. 

n 


(125) 


■/  / 

y  X 
■'o  o 


'’n'V'o’ 


II 


Q“  (k  ,k  )e 
EJyy  x’  y 


jO<x(x-Xo)+ky(y-yo) ) 


dk  dk  dx  dy 
X  y  o  0 


Substituting  (125)  into  (39)  yields 


.  ~r.  -■ 


j(kx(x-Xo)+ky(y-yo)) 


y  ^  ’'o 


dk  dk  d  d  dxdy 
X  y  X  y 
o  0 


=  -//  Qn,  (k  .K  )J^  (X  .X  )J^(k  ,k  )dk  dk 
mn  EJyy  x’  y  m  x’  y  n  x’  y  x  y 


(127) 


where 


r  r  -j(k  x  +k  y  ) 

J^(k  ,k  )  =  /  /  J^(x,y)e  ^  °  ^  °  dx  dy 

n  x  y  ■'y''xn  ''  o-^o 

(k  ,k  )  =•  j  J  J^(x,y)e  ^  ^  dxdy 

m  X  y  ■*  y'  X  m  *•' 


(128a) 


(128b) 


Thus,  the  six  fold  integration  in  (126)  can  be  reduced  to  a  double 

integral  by  Fourier  transforming  the  basis  and  test  currents 

analytically  as  defined  in  (128).  Note  that  (128b)  is  only  valid  for 

real  currents.  The  integration  in  (127)  must  be  carried  out  numerically 

so  it  is  advantageous  to  convert  the  double  infinite  integral  in  the 

k  ,  k  plane  to  a  single  infinite  integral  and  a  finite  integral  by 
x  y 

converting  to  polar  coordinates  C6J.  Letting 
k  »  6  cosa 

X 

k  «  3  si  no 

y 

(127)  becomes 


r2K  y*  “  •  _ 

■  -L  L  Qr,  (X  ,k  )J^  (k  ,k  )J®(k  ,k  )  BdBda 
mn  ■'0  "^0  EJyy  x'  y  m  x  y  n  x  y 


(129) 
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A  further  simplification  can  be  realized  by  noting  that  the  Q  functions 


are  even  in  and  allowing  one  to  integrate  on  a  from  0  to  ir/Z  only 


[6]  as  shown  below. 


C7r/2|-®  -  »  .  _  »  - 

4-  -“'o  Jo  «EJyy<''x''<y'"»‘4'<x>''n''‘x>"'®'4‘<y>''n<V’®‘’'‘’“ 


(130) 


where 


J^(k  ,k  )  =  J^(k  )J^{k  ) 
n  X  y  n  X  n  y 


(131) 


Defining  Fourier  transform  pairs  of  the  other  currents  as 


J^,(x,y) 


Jg(x,y) 


M^^(x,y) 


J^(x,y) 


(132) 

X^^^V 

(133) 

’(k^)M^P(ky) 

(13^) 

(135) 

the  remaining  matrix  elements  can  be  written  as 


r-n/2  r®  Z  *  -  Z  *  - 

-“Jo  Jo  4yy''<x’'<y>'"l4'<x>-’^‘‘x>"'''4^><'''y>'»''«'’“ 


(137) 
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(139) 
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(140) 
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Numerical  Considerations 

The  integrals  in  (130)  and  (137)-(145)  must  be  evaluated 
numerically.  The  numerical  integration  routine  used  was  a  modification 
of  a  routine  written  by  D.  M.  Pozar  to  analyze  a  rectangular  microstrip 
antenna  fed  by  an  ideal  Impulse  current  source  [6].  This  routine 
employs  a  10  point  Gaussian  quadrature  integration  scheme.  A  special 
consideration  common  to  all  of  the  integrations  is  the  presence  of 
simple  poles  determi.'.ed  by  the  zeroes  of  (96)  and  (97)  which  are  a 
function  of  B  only.  (In  region  a  the  poles  are  determined  from  (96)  and 
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b  3 

(97)  where  d  is  replaced  by  d  and  e  by  e  . )  These  poles  correspond 
D  3  p  r* 

to  the  TE  and  TM  surface  wave  modes  respectively.  The  poles  are  first 

located  by  the  Newton-Rhapson  method,  where  the  real  part  of  the  pole 

1  /2 

always  lies  between  k  and  (e  )  k  .  The  integration  on  B  is  then 

o  r  o 

carried  out  by  numerically  integrating  on  6  up  to  within  A  =  .001  k^  of 

the  real  part  of  the  pole  (8^).  The  integration  from  3^-A  to  B^  +  A  is 

performed  analytically.  The  details  of  this  integration  are  given  in 

[6].  This  cycle  is  continued  until  the  last  pole  is  encountered.  The 

remaining  portion  of  the  integration  is  calculated  numerically  until 

convergence  is  achieved.  Convergence  is  attained  when  at  least  half  the 

elements  in  a  given  vector  change  by  less  than  .5it  over  an  integration 

interval  of  50  k  .  This  usually  occurs  somewhere  in  the  interval  100  k 
0  o 

<  3  <  200  k  . 

—  —  0 

By  the  symmetry  properties  of  Green’s  functions  and  Toeplitz 

symmetries,  which  occur  due  to  the  infinite  extent  of  the  grounded 

dielectric  slab,  only  a  single  row  or  column  of  any  matrix  must  be 

calculated.  The  above  integration  routine  was  modified  to  calculate  an 

entire  row  or  column  with  a  single  call  of  the  integration  routine. 

This  greatly  improves  the  computation  time  of  the  analysis  because  an 

entire  vector  of  about  10  elements  or  less  can  be  calculated  almost  as 

fast  as  a  single  element.  This  is  due  to  the  fact  that  all  the 

integrands  can  be  factored  into  two  terms.  One  term,  T  ^(k  ,k  ),  is 

cnst  X  y 

common  to  all  the  elements  in  the  vector  and  represents  a  vast  majority 

of  the  integrand.  The  other  term,  T  (k  ,k  ,m,n),  differs  from  element 
°  var  x’  y’  ’ 


to  element  but  is  generally  simple  in  form  and  in  most  cases  is  a  single 


cosine  term  of  the  form,  cos  (where  and  y^  are  the  y 

positions  of  the  basis  and  test  modes).  Therefore,  most  of  the 

integrand  can  be  evaluated  Just  once  at  each  point,  (k  ,k  ),  and  the 

X  o  y  o 

value  of  the  total  integrand  for  each  element  in  the  vector  at  (k  ,k  ) 

°  xo  yo 

is  then  given  by 


I.n(k  ,k  ,m,n) 
T  xo  yo 


T  ^(k  ,k  )T  (k  ,k  ,m,n) 
cnst  xo  yo  var  xo  yo 


(1^6) 


To  evaluate  T  (k  ,k  ,m,n)  N  times  (N~<10)  takes  far  less  time  than 

V  wi  1  A  \  J  Jr  \J 

to  evaluate  T^  ,  (k  ,k  )  Just  once  so  the  entire  vector  can  be 
cnst  xo  yo 

calculated  nearly  as  fast  as  a  single  element. 

Due  to  the  complexity  of  the  analysis  various  tests  were  performed 

to  assess  the  validity  of  the  numerical  results.  These  tests  included 

comparisons  between  numerical  and  analytic  results  for  various  limiting 

cases  for  which  analytic  solutions  were  possible  and  comparisons  of 

numerical  results  with  numerical  results  provided  by  D.  M.  Pozar  and  R. 

W.  Jackson  from  programs  which  have  been  proven  to  be  reliable.  Tne 

latter  type  of  comparisons  confirmed  the  accuracy  of  the  ,  z^^,  z^, 

fnn  K  X  rn 

s  3 

z^  elements.  The  y  and  y  terms  were  validated  by  comparing  the 
results  to  analytic  solutions  for  the  self  impedance  of  a  narrow  slot 
with  a  PWS  E-field  aperture  distribution  in  a  ground  plane  in  free  space 
and  in  a  very  thick  grounded  dielectric  slab.  In  the  latter  case  the 
analytic  solution  is  obtained  by  assuming  the  slot  is  in  a  homogeneous 


medium  with  an  effective  dielectric  constant,  e  =(€  +1)/2.  In  both  of 

re  r 

these  cases  the  analytic  and  numerical  results  were  in  close  agreement. 

The  c  element  was  checked  by  calculating  in  closed  form  the  coupling 
n 

between  a  short,  narrow  slot  in  a  ground  plane  with  a  cosine  aperture 
distribution  and  a  short  PWS  dipole  in  the  far  field  of  the  slot.  Again 
the  two  solutions  compared  favorably. 

Lastly,  studies  were  conducted  to  determine  the  number  and  length 
of  the  current  expansion  modes  necessary  to  obtain  stable  solutions. 
From  these  studies  it  was  determined  that  five  PWS  modes  on  the  patch,  a 
traveling  wave  mode  length  of  five  feed  line  wavelengths,  i.e.,  5Xj.,  and 
a  density  of  approximately  seven  PWS  modes  per  X  ^2  starting  X^./^  prior 
to  the  aperture  center  were  adequate  for  convergence  in  the  range  of  the 
material  parameters  studied.  In  the  case  of  the  terminated  line 
analysis  the  PWS  mode  density  was  held  approximately  constant  for  all 
stub  lengths  and  in  the  case  of  the  S-parameter  analysis  seven  PWS 
expansion  modes  are  used  on  the  feedline  from  y  =  -X^./^  to  y  =  X^/4. 


CHAPTER  III 


RESULTS 


In  all  subsequent  input  impedance  plots  certain  conventions  are 


followed.  A  measured  impedance  locus  is  denoted  by  a  solid  curve  and 


the  frequency,  given  in  megahertz,  is  marked  on  the  locus  in  25  MHz 


intervals.  Calculated  input  impedance  points  are  given  in  25  MHz  steps 


and  are  connected  by  a  best  fit  dashed  circle.  Both  the  measured  and 


calculated  input  impedances  are  phase  referenced  to  the  center  of  the 


aperture.  Empirical  phase  referencing  was  achieved  by  constructing  a 


reference  microstrip  line  of  the  appropriate  length  and  shorting  the  end 


to  the  ground  plane  with  a  wide  strip  of  copper  foil. 


Figure  3.1  shows  a  measured  input  impedance  plot  and  various 


calculated  points.  The  antenna  and  feed  dielectric  substrates  were  both 


Oak  601  with  a  nominal  dielectric  constant  of  2.5*^  and  a  measured 


thickness  of  .16  cm.  The  antenna  dimensions  were  ^.0  cm  (resonant 


dimension)  by  3.0  cm.  The  slot  was  1.12  cm  long  by  .155  cm  wide.  The 


microstrip  feed  line  was  .4112  cm  wide  and  extended  2.0  cm  beyond  the 


center  of  the  aperture.  Three  sets  of  calculated  data  points  are  shown. 


The  points  defined  by  open  squares  were  obtained  assuming  a  cosine 


aperture  E-field  distribution  in  eqn.  12)  whereas  the  open 

ap  3  P 


and  filled  circles  result  from  a  PWS  aperture  E-field  distribution.  The 


input  impedance  data  indicated  by  open  circles  was  calculated  using  the 


r"-.  •.*,  ■-■V-H'.v.' 


Figure  3.1  Measured  versus  calculated  input  impedance  of  an  aperture 
coupled  microstrip  antenna  with  low  dielectric  constant  feed  sub¬ 
strate.  □  »  cosine  E-field  aperture  distribution,  direct  method  of 
calculation,  O  =  PWS  E-field  aperture  distribution,  direct  method 
of  calculation,  •  =  PWS  g-field  aperture  distribution,  S-parameter 
method  of  calculation.  (e°*2.54,  dj^=.16cm,  1^*4. Ocn.,  Wp*3.0cm, 

X  »0.0cm,  y  *0.0cm,  L  =l.l2cm,  W  ». 155cm,  £  =2.54,  d  *.16cm, 
os  •'os  ap  ap  r  a 

Wp*. 442cm,  L  =2. 0cm) 


E 


first  scheme  discussed  in  Chapter  II  where  the  feed  line  is  terminated 
in  an  open  circuit  (to  be  called  the  direct  method).  The  filled  circles 
were  calculated  by  the  S-parameter  method  discussed  in  the  second 
section  of  chapter  two.  In  the  case  of  the  PWS  aperture  distribution 
the  parameter  in  eqn.  12  was  chosen  to  be  the  arithmetic  average  of 
the  effective  wave  numbers  for  the  feed  and  antenna  regions,  i.e.,  k 

ap 

(k^  k^)/2.  This  results  in  an  aperture  E-field  distribution  which  is 

e  e 

roughly  triangular  for  the  electrically  short  apertures  used  in  these 
antennas. 

The  PWS  aperture  distribution  yields  significantly  closer  agreement  with 
measured  values  than  the  cosine  distribution.  The  results  obtained  with 
the  direct  method  and  the  S-parameter  method  are  approximately 
equivalent  but  the  latter  appears  to  match  the  measured  results  slightly 
better.  Although  the  comparison  between  analytic  and  experimental 
results  is  good,  even  better  agreement  could  be  obtained  by  varying  the 
dielectric  constant  within  the  manufacturer's  specifications  and 
adjusting  slightly  the  reference  plane  which  is  empirically  only 
approximately  known. 

To  further  confirm  the  validity  of  the  PWS  aperture  distribution 

the  same  comparison  as  above  is  shown  in  Figure  3-2  for  an  antenna  with 

a  high  dielectric  constant  feed  substrate.  The  antenna  substrate  and 

dimensions  are  the  same  as  above  but  buroid  6010.2  (e  «10.2,  d  -.1325cm) 

r*  d 

was  used  for  the  feed  substrate.  The  aperture  dimensions  in  this  case 
were  .913  x  .108  cm  and  the  feed  width  and  stub  length  beyond  the 


Figure  3.2  Measured  versus  calculated  input  impedance  of  an  aperture 
coupled  microstrip  antenna  with  high  dielectric  constant  feed  sub¬ 
strate.  a*  cosine  E-field  aperture  distribution,  direct  method 
of  calculation,  O  =  PWS  E-field  aperture  distribution,  direct  method 
of  calculation,  •  =  PWS  E-field  aperture  distribution,  S-parameter 
method  of  calculation. 

(e^=2.54,  dj^=.16cm,  Lp=4.0cm,  Wp=3.0cm,  x^^=0.0cm,  y^g=0.0cm, 

L  *. 913cm,  W  =.  108cm,  £^=10.2,  d  =.1325cm,  Wp=. 116cm,  L  =1. 107cm) 
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plotted  on  an  admittance  chart.  This  is  to  be  expected  based  on  the 
proposed  model.  Figure  3.6  shows  a  circuit  model  for  the  antenna  of 
Figure  3.1.  From  the  S-parameter  analysis  the  reflection  coefficient  at 
the  end  of  open  ended  microstrip  line  with  the  dimensions  given  in 
Figure  3.1  at  2.2  GHz  equals  .99-j.12  or  equivalently  Z  ^  (bar  indicates 
normalized  impedance)  equals  0.6-jl6.2.  The  stub  length  beyond  the 
center  of  the  aperture  at  2.2  GHz  equals  .2l4Xj.,  where  X  is  the  feed 
line  wavelength  ( X  j.=2Tr/k^) .  Therefore,  the  normalized  impedance  looking 
into  the  stub  at  the  center  of  the  aperture,  is  .002-j.17  or 
Zgtub  ~  Presumably,  the  capacitive  reactance  cancels  the  series 
inductive  reactance  over  the  bandwidth  of  the  parallel  RLC  network 
yielding  the  bottom  circuit  of  Figure  3.6.  The  input  admittance  of  this 
circuit  follows  a  constant  conductance  circle. 

Three  parameter  studies  were  carried  out  both  emperically  and 
analytically  with  the  antenna  of  Figure  3.1.  These  parameters  were 
selected  due  to  the  reasonable  ease  and  accuracy  with  which  they  could 
be  varied  experimentally. 

The  effect  of  stub  length,  i.e.,  the  length  of  feed  line  beyond  the 
center  of  the  aperture,  on  the  input  impedance  locus  was  examined.  The 
measured  loci  versus  stub  length  are  shown  in  Figure  3.7.  As  the  length 
of  the  open  circuited  stub  is  decreased,  from  an  initial  dimension  of 
.214  X^,  the  input  impedance  at  any  given  frequency  moves 
counterclockwise  along  a  constant  resistance  circle  towards  the  open 
circuit  point  on  the  Smith  Chart.  Since  the  aperture  looks  like  a 


Figure  3.8  Measured  versus  calculated  input  impedance  loci  as  a 
function  of  stub  length.  (See  figure  3.1  for  antenna  parameters.) 

(a)  Calculated  points  obtained  with  direct  method. 

(b)  Calculated  points  obtained  with  S-parameter  method. 


Figure  3.8  continued 
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series  impedance,  only  the  imaginary  part  of  the  input  impedance  at  any 
frequency  changes  as  the  length  of  the  open  circuited  stub  varies. 

A  comparison  of  the  above  result  with  calculated  loci  is  shown  in 
Figure  3.8-  Only  every  other  curve  of  Figure  3.7  is  shown  in  order  to 
improve  the  clarity  of  the  comparison.  Input  impedance  data  calculated 
by  the  direct  method  is  displayed  in  Figure  3.8a,  whereas  the  S- 
parameter  method  was  used  to  produce  the  analytic  impedance  plots  in 
Figure  3.8b.  In  both  cases  the  analytic  curves  approximate  the  measured 
loci  reasonably  well. 

Another  empirical  study  conducted  with  the  antenna  of  Figure  3.1 
was  to  vary  the  patch  position  relative  to  the  aperture.  Measured  and 
calculated  plots  are  given  in  Figure  3-9  correspond xng  to  movement  of 

A 

the  patch  in  the  y-direction,  i.e.,  along  the  resonant  dimension  (see 
Figure  2.2a).  Calculated  points  utilizing  the  direct  method  and  the  S- 
parameter  method  are  plotted  in  Figure  3.9a  and  3.9b  respectively.  The 
agreement  is  good  in  both  cases.  The  coupling  factor,  as  defined  by  the 
radius  of  the  resonance  circle,  is  greatest  when  the  patch  is  centered 
over  the  aperture  and  drops  significantly  as  the  patch  is  moved  in  the 
y-direction.  This  is  in  accordance  with  Pozar's  [1]  simple  model  for 
this  antenna  based  on  Bet he  hole  theory  and  the  cavity  model.  In 
addition,  as  the  patch  is  offset  in  the  y-direction  the  centers  of  the 
resonant  loops  move  approximately  in  a  straight  line  towards  the  edge  of 
the  Smith  chart  just  to  the  inductive  side  of  the  short  position.  This 
is  probably  due  to  the  fact  that  when  the  patch  is  offset  by  a  large 


Figure  3.9  Measured  versus  calculated  input  impedance  loci  as  a 
function  of  patch  offset  in  the  direction  of  resonance.  The  schematic 
above  shows  the  relative  position  of  the  patch  to  the  slot  in  each 
case.  (See  figure  3.1  for  other  antenna  parameters.) 

(a)  Calculated  points  obtained  with  direct  method. 

(b)  Calculated  points  obtained  with  S-parameter  method. 


amount  the  structure  looks  like  a  small  aperture  in  a  ground  plane  which 
is  inductive. 

In  contrast  to  offsetting  the  patch  in  the  y-direction,  lateral 
movement  of  the  patch  in  the  x-direction  causes  little  change  in  the 
coupling  factor  provided  the  entire  slot  remains  under  the  patch.  From 
the  measured  data  in  Figure  3.10  it  can  be  seen  that  the  coupling  factor 
actually  increases  as  the  edge  of  the  patch  aligns  with  the  edge  of  the 
slot  and  then  monotonically  decreases  as  the  slot  emerges  from  under  the 
patch.  Figure  3.11  shows  calculated  loci  versus  x-offset  of  the  patch. 
Initially  the  coupling  factor  remains  constant,  however  at  the  point 
where  the  coupling  factor  increases  empirically  it  decreases 
analytically.  This  causes  a  significant  discrepancy  between  measured 
and  calculated  results  for  all  large  patch  offsets  in  the  x-direction. 
This  disagreement  is  not  unexpected  since  our  model  utilizes  only  one 
mode  in  the  aperture.  A  single  aperture  mode  makes  the  analysis 
numerically  more  tractable  but  cannot  account  for  skewing  of  the 
aperture  electric  field  distribution  as  the  patch  is  offset  in  a 
direction  parallel  to  the  long  dimension  of  the  slot.  In  addition,  the 
patch  current  is  assumed  uniform  in  the  x-direction  which  may  not  be 
adequate  for  large  offsets  in  that  di'^ection. 

All  subsequent  parameter  studies  are  based  solely  on  analytic 
results.  The  direct  method  of  calculating  input  impedance  was  selected 
rather  than  the  3-parameter  method  for  these  studies  because;  (1)  the 
agreement  .between  the  two  methods  is  within  the  limits  of  measurement 


Figure  3.11  Calculated  (direct  method)  input  impedance  loci  as  a 
function  of  patch  offset  in  the  direction  orthogonal  to  resonance. 
The  schematic  shows  the  relative  position  of  the  patch  to  the  slot 
in  each  case.  (See  figure  3.1  for  other  antenna  parameters.) 


accuracy,  (2)  the  direct  method  is  slightly  faster  except  in  the  case  of 
stub  length  variation  studies,  (3)  all  currents  for  the  tuned  antenna, 


i.e.,  feed  line  terminated  in  an  open  circuit,  are  obtained  explicitly 
with  the  direct  method  as  a  function  of  the  parameter  varied. 

The  long  dimension  of  the  aperture  was  varied  to  obtain  the  curves 
given  in  Figure  3.12.  The  antenna  dimensions  are  given  in  the  Figure 
legend  and  are  very  similar  to  the  dimensions  of  the  antenna  of  Figure 
3.1.  As  the  aperture  length  is  reduced  the  radius  of  the  impedance 
circle  decreases  and  the  center  of  the  circle  moves  towards  the  short 
circuit  location.  The  resonant  frequency  versus  aperture  length  is 
plotted  in  Figure  3.13.  The  resonant  frequency,  which  in  this  case  is 
also  the  minimum  VSWR  frequency,  decreases  with  increasing  slot  length. 
Also  plotted  in  Figure  3.13  is  the  input  impedance  at  resonance  versus 
slot  length  which  can  be  used  to  approximately  determine  the  slot  length 
required  to  achieve  critical  coupling  and  the  correspond  ing  resonance 
frequency.  In  this  case  the  aperture  length  which  yields  critical 
coupling  is  1.09  cm  at  a  resonant  frequency  of  2.233  OHz.  For 
comparison  the  resonant  frequency  of  this  antenna  based  on  the  cavity 
model  is  2.306  GHz  [3]. 

It  is  also  of  interest  to  examine  the  influence  of  feed  substrate 
dielectric  constant  and  thickness  on  the  input  impedance,  since  the  feed 
substrate  will  be  electrically  thick  in  the  intended  application  of  the 
antenna.  As  dielectric  constant  and  thickness  are  varied  in  these 
studies  the  feed  line  width  and  stub  length  are  modified  to  maintain  a 


Figure  3.12  Calculated  (direct  method)  input  impedance  loci  as  a 
function  of  aperture  length  (long  dimension).  The  other  antenna 
parameters  are: 

^-“2.54,  d,=.16cra,  L  =4.0cm,  W  =3. 0cm,  x  =0.0cm,  y  =0.0cm, 

‘  °  a  "  P  os 

W  =  .11cm,  e  =2.54,  d  =.16cm,  W^=. 495cm,  L  =2. 0cm 


Slot  LongtM  (cb) 


Figure  3.13  Resonant  frequency  and  input  resistance  at  resonance 
versus  slot  length  (data  from  figure  3.12). 
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characteristic  impedance  of  50  2  and  a  stub  length  of  .22  \ All  other 
antenna  parameters  were  held  constant  and  are  given  in  the  figure 
legends. 

The  dielectric  constant  variation  is  shown  in  Figure  The  key 
features  are  the  increase  in  the  coupling  factor  and  the  invariance  of 
the  resonance  frequency  with  increasing  dielectric  constant.  The 
increase  in  the  coupling  factor  is  probably  due  to  the  fact  that  the 
electric  length  of  the  slot  is  increasing  as  the  dielectric  constant  of 
the  feed  increases. 

The  last  parameter  study  performed  was  to  increase  the  thickness  of 
the  feed  substrate  of  the  antenna  of  Figure  3-1**  in  the  case  of  e^=10.2. 
As  the  distance  between  the  feedline  and  aperture  increases,  the 
coupling  factor  decreases  as  can  be  seen  in  Figure  3.15.  As  with  the 
dielectric  constant  variations,  the  resonant  frequency  is  unchanged  with 
changes  in  substrate  thickness  over  the  range  studied. 
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Figure  3.15  Calculated  (direct  method)  input  impedance  loci  as  a 
function  of  feed  substrate  thickness.  The  tabular  data  above  gives 
the  feed  line  width  and  stub  length  used  in  the  analysis  to  main¬ 
tain  a  50Q  characteristic  impedance  and  a  stub  length  of  .22X, 
for  each  value  of  d  .  The  other  antenna  parameters  are: 
b  ^ 

e_=2.54,  d.=.16cm,  L  =4. 0cm,  W  =3. 0cm,  x  =0.0cm,  y  =0.0cm, 
r  b  ’  p  ’  p  os  ■'os  ’ 

L  -l.Ocrn,  W  =.llcm,  e^=10.2. 
ap  ap  r 


CHAPTER  LV 


CONCLUSION 

A  rectangular  microstrip  antenna  coupled  to  a  microstrip  line 
through  a  small  rectangular  aperture  in  the  ground  plane  has  been 
analyzed  by  the  moment  method.  In  the  case  of  a  feed  line  terminated  in 
an  open  circuit  the  input  impedance  is  determined  from  the  amplitude 
coefficient  of  a  reflected  traveling  wave  current  mode  on  the  feed  line. 
In  the  case  of  an  infinite  feed  line  the  S-parameters  are  calculated  by 
including  both  reflected  and  traveling  wave  current  modes  on  the  feed 
line.  These  amplitude  coefficients  are  obtained  directly  from  the 
moment  method  solution. 

The  analysis  has  been  verified  by  comparison  with  measured  input 

impedance  data  for  an  antenna  with  a  low  dielectric  constant  (e  =2.5^) 

r 

feed  substrate  and  one  with  a  high  dielectric  constant  (e^=10.2)  feed 
substrate.  In  the  former  case  empirical  parameter  studies,  which 
involved  varying  the  length  of  the  feed  line  beyond  the  aperture  and 
lateral  displacements  of  the  antenna  with  respect  to  the  aperture,  were 
carried  out  and  compared  with  anlytical  studies  to  further  validate  the 
theory.  With  the  exception  of  patch  offsets  orthogonal  to  the  resonant 
dimension  of  the  patch  the  analytical  and  empirical  results  were  in  good 
agreement.  From  the  above  data  and  S-parameter  results  an  equivalent 
two-port  lumped  element  circuit  is  proposed  for  the  aperture  backed  by  a 


microstrip  antenna  near  resonance.  The  circuit  is  a  series  impedance 


consisting  of  a  parallel  RLC  circuit  in  series  with  an  inductor. 
Finally,  analytical  parameter  studies  have  been  performed  to  determine 
the  effect  of  aperture  length  and  feed  substrate  dielectric  constant  and 
thickness  on  input  impedance  and  resonant  frequency. 

Further  development  of  the  proposed  equivalent  circuit  for  this 
antenna  should  be  undertaken  in  the  future.  The  circuit  elements  could 
be  found  by  determining  the  stub  length  that  yields  an  admittance  locus 
which  follows  a  constant  conductance  circle.  The  series  reactance  due 
to  L'  would  then  be  the  negative  of  the  reactance  looking  into  the  open 
circuited  stub,  from  which  L'  could  be  determined.  The  resistance,  R, 
of  the  model  would  simply  be  the  inverse  of  the  input  conductance 
(constant  under  the  condition  stated  above)  near  resonance  and  L  and  C 
could  be  determined  from  the  input  susceptance  at  two  frequencies  around 
resonance.  By  performing  parameter  studies  and  examining  the  effect  on 
the  equivalent  circuit  element  values,  the  antenna  parameters  which  most 
strongly  influence  each  circuit  element  might  be  elucidated.  Another 
area  for  future  worK  is  to  determine  the  far-field  antenna  patterns  from 
the  currents  on  the  patch,  aperture  and  feed  line.  Of  particular 
interest  would  be  the  relative  magnitude  of  the  lobe  on  the  feed  side 
versus  that  on  the  antenna  side.  Finally,  improved  programmming 
techniques  should  be  sought  to  reduce  the  CPU  time  required  to  carry  out 
the  analysis.  Most  importantly,  this  would  make  it  feasible  to  expand 
the  aperture  magnetic  current  in  more  than  one  basis  mode. 
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